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and discuss in detail the effective coupling constants for the marginally long ranged 
dipole-dipole interaction. 
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1. Introduction 



The achievement of Bose-Einstein condensation in dilute ultracold gases JTJ also renewed 
the theoretical interest in fermionic many body systems |2], [| [|, e.g. their superfluid 
properties |j[ 0, |, |TT| . Experimental successes in obtaining degeneracy in Fermi 
vapours |TT|, [13], |15[] further stimulated the interest in confined ultracold Fermi 
gases. Recently, successful production of large samples of fermionic lithium (10 6 - 10 7 
atoms) with long lifetimes ( 10 2 sec.) was reported fl6|, [17| . 

The present situation in this field is characterized by intensive search for a verification 
of the presence of a superfluid BCS-like phase. Strong attractive interaction between 
the fermions is needed for a transition temperature not too far below the temperature 
of Fermi degeneracy. It is not only the absence of s-wave scattering between identical 
fermions which causes problems, even the intrinsic interaction between non-identical 
fermions is generally too weak. Strongly interacting degenerate Fermi gases now can be 
produced either by mixing bosons with fermions [pl| [19|] or by magnetically inducing a 
Feshbach resonance in a two-species mixture of fermions [fj(], [H], [2^, [23[| . 
Here we study the interaction of both identical and non-identical fermions in a 
one-dimensional harmonic trap where no superfluid transition can occur at finite 
temperatures. The interactions considered in the present paper are limited to scattering 
in a single channel when energy is exchanged between the scattering particles and the 
harmonic trap. Moreover, a large interaction energy would violate the criterion for a 
one-dimensional regime as is discussed in Section 2. 

Using microtrap technology |f24| , [25| , |2~6 , [27], |28], it will become possible in the near future 
to produce a neutral ultracold quasi one-dimensional Fermi gas in a highly elongated 
trap. 

The confinement of a trapped ultracold gas can be realized by a harmonic external 
potential. Theoretical work on interacting fermions in a one-dimensional harmonic trap 
has been published in |30], |31 |. It is based on concepts of Luttinger liquid theory |32| (for 
reviews see [[33], |35|]). Detailed knowledge of the coupling constants will be necessary 
for understanding the ground state properties of this novel confined one-dimensional 
Fermi gas. The coupling constants will decide which type of fluctuations are dominant, 
density wave, component wave, or pairing fluctuations. 

Solvable models require specific forms of the interaction coefficients, especially the 
exclusion of momentum non-conserving collision events in the harmonic trap (which 
evidently are present here). We study the coefficients of the pair interaction in order to 
find out to what extent physical interactions have at least approximately the behaviour 
required by the exactly solvable models. 

We construct an effective one-dimensional potential from a physical interaction potential 
using the simplest conceivable realization of the projection from three dimensions to 
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one dimension. We set up the effective interaction operator by taking into account 
only interaction coefficients of the transverse ground state and study the effects of the 
one-dimensional geometry on the form of the interaction coefficients. 
In many cases, identical spin polarized fermions experience only a weak residual 
interaction because s-wave scattering is forbidden. This restriction does not hold for a 
two- component system of spin polarized fermions where significant interactions between 
the components are possible. For instance, the dipole-dipole interaction [SB] can become 
relevant, especially in the case of polar molecules |S7| ]. 

Postulating the existence of an effective local one-dimensional potential with a Fourier 
transform, we demonstrate that three main scattering processes dominate the interaction 
coefficients in one dimension. They are identified as two forms of forward and one 
form of backward scattering. They determine the coupling constants in theories of the 
interacting one-dimensional degenerate Fermi gas. 



Recently, a two-component Fermi gas has been realized |38| with different hyperfine 
states of 40 K. For simplicity, however, we assume the same longitudinal trapping 
frequency for both components. 

The paper is organized as follows. Section 2 develops the theory for the interaction 
coefficients. It includes a discussion how the one-dimensional regime may be reached. 
In Section 3 we evaluate the interaction coefficients for an arbitrary interaction of 
finite range in terms of the zeroth and second moments of the effective one-dimensional 
interaction potential. In Section 4 we present the results for the marginally long-ranged 
dipole-dipole interaction. In both cases we assess the relevance of our findings for one- 
and two-component Fermi vapours in one dimension. 

2. Theory 

Starting from the case of identical atoms in the same hyperfine state, the pair interaction 
operator V3 for a many particle system in a trap is 

^3 = V(m,p 1 q,n)cj ri c q cp~c n , (1) 

with interaction coefficients 

V(m,p;q,n) = J d 3 xdV^(x)^(x')l/ 3 (x-x')^(x)^(x'). (2) 

V5j(x — x') denotes the two-particle interaction potential while the quantum numbers 
m x ,m y ,m z etc., abbreviated by m, characterize the single particle bound state in the 
potential of the trap. 
In our case, 

1 1 

V(x) = -m A u 2 z 2 + -m A u 2 ± (x 2 + y 2 ) (3) 
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is the potential of the highly elongated (trap frequencies uj± ^> ue) axially symmetric 
harmonic trap. The atomic mass is denoted by m^. The non-interacting system has 
single particle levels 

huj n = hu e (n z + ^ + %u±_{n x + n y + l), n x , n y , n z = 0,l,... (4) 

and wave functions which are products of the wave functions of the one-dimensional 
harmonic oscillator, e.g., for the coordinate z in the elongated direction of the trap 

(n = n z ) 

M*) = (^) exp(-aV/2) Kn(az). (5) 

H n denotes a Hermite polynomial. Intrinsic length scales of the trap are i = a^ 1 
with a 2 = rriAUJi/h and £± = a J 1 with a\ = m^uj^/fi^ the spatial extensions of the 
longitudinal and transverse ground states, respectively. 

We consider a system of 1 <C iV < uj^/uj^ fermionic atoms. At T = 0, they fill the lowest 
N states well below the first transverse excited state. Note that each single particle level 
is occupied by at most one spin-polarized fermion. Then the unperturbed Fermi energy 
is 

e F = fiw e ( N -\)- (6) 

The transverse ground state energy hu± is taken from now on as zero of energy. 
Neglect of all terms with transverse quantum numbers not equal to zero defines our 
quasi one-dimensional model. The validity of this approximation is discussed below. 
From now on we write m for m = (0,0, m z ) etc. and use an effective one-dimensional 
potential Vi e fr which results from averaging the physical potential V3 over the transverse 
ground state 

V leS (z - z') = [ d 2 xd 2 3 ;> ±0 (x) 2 ^o(y)V±o(a:') 2 ^o(2/ / )V3(x - x'). (7) 



For example, a three-dimensional contact potential V^(x) = g8^(x) (which contributes 
only for bosons or non-identical fermions to V3) gives 

V leS (z) = g^S(z). (8) 

Instead of computing V(m,p] q,n) from the real space potential equation (|7|) we use in 
our calculations below a Fourier version of equation (0) which reads 

/ /, , ,n / n ( minfm, q)\ minfra, p)\ \ ^ 2 . , 

V{m,p; q, n) = cos(vr(|g - m\ - \n - p\)/2) 2 )— (9) 



max(m, q) \ max(n,p)! 
X ^L [™ dvv^- m ^ n - p ^ 2 e~ v V lcfi (k 2 = 2a 2 v)L {lq r^\(v)L (ln ~ pl \ [ 

2^j- Jq \ II ' tmn(q,m) \ 1 mln(?^,p) , 
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provided the Fourier transform Vjeff(fyl) of Vx e s(z) 
V leS {k\\) = J dzexp(ik\\z)V le{i (z) 



exists. The z-component of the Fourier variable k is denoted by k\\. denotes a 
Laguerre polynomial. 

We generalize to the case of two components and discuss the conditions under which the 
one-dimensional model defined above is valid. First of all we exclude thermal excitations 
of the transverse levels by demanding ksT <C Huj±. Consider now 2N fermionic atoms 
in an elongated harmonic trap which are distributed equally into 2 hyperfine states. 
Atoms of different components interact through a three-dimensional contact potential 
g5^ 3 \x. — x'). Interaction between atoms of the same component is much weaker and 
is therefore neglected here. The relation iV < ou±/uj£ guarantees that all the atoms can 
be accomodated in the non-interacting ground state. We use the mean field estimate 
ng for the interaction energy per atom denoting the density in three dimensions of each 
hyperfine component by n. Following D. Gangardt and and G. Shlyapnikov this 



energy must be much smaller than the transverse excitation energy hu± to be in the 
one-dimensional regime. Hereby it is assumed that the width of the transverse ground 
state of the trap greatly exceeds the three-dimensional scattering length a = mg/Anh 2 
for the scattering of fermions in different hyperfine states |40|| . N atoms in the transverse 



ground state of the trap occupy a volume of ~ 2tc a 1 a 2 \/2N and the above criterion 
may be written as 

aaV2N < 1. (10) 

Note that the aspect ratio (uje/uj±) 1 ^ 2 does not enter the criterion ([10| ) explicitely. It is 
not only a sufficiently large ratio of LU±/u}g but the reduction of iV which leads into the 
one-dimensional regime provided the Fermi level is always below hu±. 
As an example we consider fermionic potassium atoms in 2 different hyperfine states in 
an elongated harmonic trap with typical parameters: Ui/u±_ = 10~ 4 , a = 10 3 cm -1 , and 
a = 8 nm. The transverse ground state of the trap be half filled by each component. 
Then the ratio between the energies of interaction and transverse excitation is as small 
as 0.05. In what follows we assume that the above criterion is always fulfilled. 
We continue to discuss the interaction coefficients V(m,p; q,n). From equation (|7|) the 
following features may be deduced: If the atoms interact by the marginally long ranged 
dipole-dipole interaction which has no intrinsic range, the effective one-dimensional 
potential has a finite range dm ~ £±. This is shown in detail in Section 4. We note 
that a similar calculation may be performed for the Coulomb interaction which has an 
infinite range and results in dm ~ 3£j_. If instead the interatomic potential has a finite 
range d^r> as e.g. the van der Waals interaction then we find dm ~ d^,. 
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In addition, the integrand in equation (9) contains further functions including Laguerre 
polynomials. This weighting provides a soft momentum cut-off at k > 2k-p or v > AN, 
provided we consider the physically relevant case m~n~p«g~iV^>l, but 
\m — n\ <C N, \m — q\ <C N, etc.. It is here, where the fermionic nature of the quantum 
gas becomes relevant. 

There are 24 coupling coefficients V(m,p; q, n) for any set of four distinct integers. Due 
to the symmetries [m «-> q), (n «-> p), and ([m, q] ^ [n,p]) only three of them can be 
different. 

It is one of the aims of the present work to study how far atoms confined to a 
one-dimensional harmonic trap and interacting by physical potentials have interaction 
coefficients V(m,p; q,n) which behave as 



Vym, p; q,n) ~ V a 5 m - qtn - p + Vb 8q-m,n-p + V c 5 m + q ^ n + p . (11) 

For such an interaction an exact solution of the model was found in previous work |30[ . 
Fermions in a box have exactly the required behaviour as is shown in the appendix. This 
is a consequence of local translational invariance inside the box. In the one- dimensional 
harmonic trap, however, the situation is different. Only in a region around the center of 
the trap, say at \z\ < 5, the atoms nearly feel no external force. Here the single particle 
states are superpositions of plane wave states exp(ik n z) with k n = ±a^2n — 1. Because 
the relevant states are near the Fermi edge, we have \k n \ ~ k-p. According to equation 
(H), incoming states {n, q} are transformed into outgoing states {p, m} in the collision 
process. Denoting a state with k n w — k? by (— n), three types of collision processes 
may be discriminated 

a-type: {n, q} -> {p,m} } 

b-type: {n, -q} -> {p,-m}, (12) 
c-type: {n,—q} — > {— p, m}. 

If the atoms move in the same direction, only forward scattering (Ap -C hkp, type a) 
may occur. If the atoms move in opposite directions, forward scattering (Ap hk-p, 
type b) as well as backward scattering (Ap w 2hkp, type c) are possible. The first two 
cases were considered in the first two papers of [3(| . The last one requires an extension 
of the bosonization method, which is the aim of the third paper [ 30(| . The couplings V a , 
Vb, and 14 are the analogues of the Luttinger model couplings g 4 , g 2 , and g±, respectively 



34, 35 



Because the quantum numbers m,p,q,n denote energies in units of hue, momentum 
conservation can only be approximately fulfilled. If e.g. in the c-type collision process 
m + q — p — n = then, linearizing with respect to m — N etc. gives 

\k m \ + \k q \-\k p \-\k n \=0(N- 1 ), (13) 
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and one would conclude that the coefficients V a ^ yC would dominate more clearly with 
increasing N. 

On the other hand, the larger N, the larger is the fraction of atoms in the regions 
\z\ > 5 where momentum is not conserved (recall that the number of atoms in the 
central region \z\ < 5 is only oc V~N) and this would enhance the momentum non- 
conserving coefficients. 

From the discussion below it will turn out that for short range interactions 
(approximated by a contact potential) coefficients of type a,b, or c dominate over the 
coefficients with indices which do not fulfill one of the three momentum conservation 
relations above. Moreover, we will show that backscattering is strongly suppressed for 
interactions with a range of more than the interparticle separation. 

3. General Discussion 

Before going to specific interactions, we address some general questions related to the 
above results. 

Interactions produced by a one-dimensional contact potential will henceforth be called 
"s-wave". Such a contact potential produces interaction coefficients, which are fully 
symmetric in all their arguments. Therefore they do not contribute to the interaction 
between identical fermions. We demonstrate how a contact contribution is subtracted 
from the interaction coefficient 

V(m,p;q,n) = J dz dz' ip m (z)^ q (z) V(z - z) ip p (z')ip n (z'). (14) 

The one-dimensional potential V(z) could, for instance, be an effective one-dimensional 
potential calculated from a physical potential via equation (0). If V(z) decays faster 
than the product i/) p (z')ip n (z') oscillates (for p « n w N the wavelengths are about 
n/k-p ~ 7r£j_) one can expand ip p (z')ip n (z') at z' — z up to second order and gets 

V(m,p;q,n) = (V) J dzip m (z)ip q (z) ip p (z)ip n (z) 
If d 2 

+ - (z 2 v) j dzij m (z)M z ) ^2 W> P ( z )M z )] + ■■■ 

= V^ 8 '(m,p; q, n) + V^ p '(m,p; q,n) + . . . (15) 

We call the first term in equation (15), which corresponds to a contact interaction, 
s-wave") and the second term ("p-wave") contribution to the interaction 
coefficient. Higher order "partial waves" can be defined analogously. 
A sufficient condition for the validity of the expansion is "C 1. Here we use the 

notation (V) = J dz V(z) and assume that the first moment (z V) of V vanishes because 
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V is an even function. The moments of V may be expressed in terms of the Fourier 
transform V(k) of the potential and its second derivative V"(k) at k = 

(V) = f dzV(z) = V(k = 0), (16) 



(z 2 V) = Jdzz 2 V(z) = -V"(k = 0), (17) 

provided the integrals converge. The Fourier variable k conjugate to z in this case 
corresponds to k\\ as explained above. 

We note that the expansion equation (15) is equivalent to inserting the approximation 
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V(z) » (V) 8{z) + I (z 2 V) ^-5(z) + . . . (IS 



2 N ' dz 



into equation (|14|) . 

In a contact interaction, the momentum transfer is uniform: Both forward scattering 
(k <C fcp) as well as backward scattering (k ~ 2kp) are equally present in the potential 
and hence contribute equally to the first term in equation (15). In "p-wave" scattering, 
however, the second derivative leads to a suppression of the forward contribution. 
In the Fourier domain, equation (|18|) reads 

V(k) = V{k = 0) + i V"(k = 0)k 2 + 0(fc 4 ), (19) 

provided V(k = 0) and V"(k = 0) exist. Again, the first term is called "s-wave" and 
the second " p-wave" . 

The dominant interaction between different components in a mixture of Fermi gases 
with short range interaction is the "s-wave" part. It is characterized by interaction 
coefficients, which are fully symmetric in its arguments. 

In this case, inserting only the "s-wave" part V(k = 0) of equation flT9| ) into equation 
(9), the v-integration can be performed exactly, leading to the result: 



V is \m,p; q,n) = a(V) 



\ 1 /2 

cos(7r(p + q — m — n)/2) ( T(m + 1) \ 



V2~n \r(p+l)T(q+l)r{n + V 

T(a 1 )T(a 2 ) 



T(a 3 )T(m-q+l) 



F 2 (—q, a4,a 2 ;a 5 ,m- q+ 1; 1) (20) 



with the abbreviations ax = (p + q + n — m + l)/2, a 2 = (m + p — q — n + l)/2, 
a 3 — ip + Q — rn — n + l)/2, a 4 = (m + n — p — q + l)/2, a 5 = (m — n — p — q + l)/2. 
The generalized hypergeometric function ^F 2 reduces here to a finite series 



The coefficients V^(m,p;q, n) have arguments m,p,q,n ~ iV while the differences of 
the arguments are much smaller than N as was explained above. In this case, the relative 
variation of the arguments of the functions which determine (m, p; q, n) according to 
equation (20) is small with respect to a variation of N but of order unity with respect 
to variations of the differences of the arguments. Numerical examination of equation 
(20) shows that V {s \ m,p; q, n) behaves indeed in this way. 

Following equation (11), we use the notation V a for V(m,p; q, n) with m + p = q + n. 
and V c are defined analogously. Specifically for the coefficients V a , Vb, and V c , we find 
from equation (20) 

K ( bc( m >P; q,n & N) oc N-°- 3 . (21) 



For the model interaction in the appendix (see equation ( A.6| ) below) we have instead 



T4 : b,c oc iV" 1 / 2 (note L pa Lp) in accordance with an argument given in |HJ. In contrast, 
the momentum non-conserving coefficients are almost independent of N. Concerning 
the neglect of all background coefficients, as in the Tomonaga-Luttinger model with 
harmonic confinement, we conclude from the present work that this model may be 
applicable to non-identical fermions if the number of atoms in the trap is limited J29 |. 



In figure 1, we show corresponding results. We have chosen the form of a histogram 
since the coefficients V(m,p; q, n) with an odd sum of arguments vanish. 
The "p-wave" interaction coefficients are obtained by inserting only the "p-wave" part 
^V"(k = 0) k 2 of equation (|T9|) into equation (9). They can also be given in closed form: 

{m,p,q,n) a [z ) ^_ ^ [r(p + l)T(q + l)T(n + 1) ) 

x n ^ K &2 ' 6s ' m ~ q + ^ ^ ^ 

T(b 3 )T(m - q + 1) 

with the abbreviations b± = (p + q + n — m — l)/2, 62 = (m+p — g — n + 3)/2, 
^3 — (p + q — Tn — n — l)/2, 64 = (m + n — p — q + 3)/2, 65 = (m — n — p — q + 3)/2. 
This formula is relevant for identical fermions with a short range interaction in a quasi 
one-dimensional harmonic trap. 

Corresponding results are shown in figure 2. Note that units are different from figure 1. 
Forward scattering is almost completely suppressed and the backward scattering 
interaction coefficients V c become dominant. The iV-dependence of ist different 
from that of V^ s \ We find 

v (p) K ^0.5 ^ (23) 

which agrees for L pa Lp completely with equation ( |A.11| ) below in the limit kpd 1 , 
even if our interaction model in the appendix is extremely crude. The dominance of 
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V c over V a and Vb persists for all N. The momentum non- conserving coefficients also 
increase with N but slightly less than V c . 

As was discussed in Section 2, the Fermi algebra removes any fully symmetric part of 
the interaction coefficients V(m,p; q,n). Inspection of figures 1 and 2 shows that not 
only the two forward scattering peaks disappear after this procedure but also almost all 
momentum non-conserving coefficients. In principle, we considered an arbitrary finite 
ranged potential V(z) here but the approximation fllBj) and flTU| ) require an interaction 
of short range to be valid. Note also that figures 3 and 4 below show that the same 
feature is present in the case of the marginally long ranged dipole-dipole interaction. 
We conclude that theories like |3(J which exactly solve the model find strongest support 
by our results if they take into account backscattering between identical fermions. 



4. Dipole-Dipole Interaction 

The real space interaction between two parallel magnetic dipoles \x is 



m!_ 3 (^m)!_^(3) (x) 



(24) 



The Fourier transform becomes 36 



(25) 



These functions are singular at x = and k = 0, respectively. 

With trivial redefinitions, all formulae can be taken over to the case of molecules 
interacting via permanent electric dipoles. 
The Fourier version of equation (0) reads 



Vleff(fc||) 



d 2 k ± 



exp(-k 2 ± /2al) V 3 (k), 



(26) 



where k», k± are the components of the wave vector k with respect to the axis of the 
trap. 

A straightforward application of equation (26) gives the effective one-dimensional dipole- 
dipole potential 



V, 



lcff DD 



(*ll 



HoH 2 a 2 ± \ 3 cos 2 $ — 1 



2tt 



K ,,2 



1 + Y e 



Ei(-^ 2 /2) 



+ sin $ 



(27) 
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with k = k\\foi.±_. Ei denotes the exponential integral and <£> is the angle between the 
dipoles and the axis of the trap. 

The effective one-dimensional potential is non-analytic at k = 0, but has the finite value 



VleSDDih = 0) = -^^(1 + COS 2 $). (28) 

For identical fermions, this value must be subtracted from equation (27), as was 
explained above. 

It is possible to calculate exactly the real space version of the effective one-dimensional 
dipole-dipole interaction: 



V lcSDD (z) = I" dh e~ ik W z Vleff dd(*||) (29) 



2tt 



8V27T 



(3 cos 2 $ 



2/na±\z\ 



2rc a± 



1 + a\z 2 ) exp (a 2 ± z 2 /2^ erfc (ai.|z|/>/2 
sin 2 $ 8{z 



For [$[ < (>) arccos 1/3, the potential is attractive (repulsive). Note that in a 
one-dimensional Fermi gas quantum fluctuations prevent the classical collapse of this 
configuration until a threshold attractive interaction strength is reached which according 



to [30] is |T4,b,c| — 0(hui). The range dm of this effective one-dimensional potential 
inside the trap follows from dma± = 1, which gives k F dm = v2 for a completely filled 
trap. The real space potential, which has a cusp at z — 0, is shown in the inset of figure 
4 for the case of dipoles oriented parallel to the axis of the trap ($ = 0). 
The large distance behaviour is 



VWd(W » al 1 ) - (3 cos 2 $ - 1) -L (30) 

This is exactly the behaviour of the initial potential equation in three dimensions (24) 
along the z-axis for two parallel dipoles. 

We now turn to the interaction coefficients for the dipole-dipole interaction. The 
analytic procedures outlined in equations (15) to ( |T9"D are not applicable to the dipole- 
dipole interaction. We thus resort to a numerical evaluation and consider the case of 
dipoles oriented parallel to the z-axis: The sequence V{m = 107, p = 103; q = 105, n) 
of interaction coefficients is displayed in figure 3, as calculated from the full equation 
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(26), i.e., for the interaction between fermionic atoms in different hyperfine states. Here 
and in figure 4 a filling factor F = Nue/u± ~ 1 is assumed. 

It is seen that the backward scattering coefficient V c is suppressed in comparison to 
V & ~ Vb (for large N, V a — > V\> hold due to the long range of the interaction. This is in 
accordance with the results for the model interaction in the box (equation (20) below, 
case 2k-pd = 0(1)). 

The iV-dependence of the coefficients V a , Vb, and V c for the full dipole-dipole interaction 
is 

K,y b oc AT ' 3 , V c ocN~ 0A (31) 

while the background coefficients which dominate more than in the case of the contact 
interaction (cf. figure 1) remain almost constant. 

In figure 4, we show the interaction coefficients for identical fermions, i.e., for the effective 
one-dimensional potential Vi e ffDD(&) — V^efr dd (^ = 0). It can be seen that in this case 
the approximation 

V{m, p; q, n) « K (n ° n_s) S m+q , n+p (32) 
is appropriate, i.e., only backward scattering is relevant. 

The iV-dependence of Vj; non ~ s ^ is different for different filling factors F = Ntue/uj±. For 
F « 1 the iV-dependence is estimated as V^ non_s " ) oc iV _a5 while for F <C 1 it becomes 

^(non-s) ^ g oth resu l ts 

agree with equation (^) as well as with equation (|A.7| ) 



(note that in the first case fcp d ~ 1 while hp d <C 1 for the second case) . The background 
coefficients grow slightly with N. 

From our results, we conclude that neglect of the background coefficients finds 
justification for the case of identical fermions in traps with a small filling factor after 
the s-wave part has been subtracted. 

Finally, we estimate V a ~ Vb for fermions in different hyperfine states, which interact 
via the magnetic dipole-dipole interaction. For m,p,q,n m N the coefficients 
V^ c (m,p;q,n) do not exceed V(N, N; N, N). With dipoles oriented parallel to the 
trap axis and using the "s-wave" approximation in equation (20) we get 

V^(N, N, N, N) « -1.34 ^ ^ ®L huJe . 
v ; 4vr h 2 y/F 



We refer to the planned highly anisotropic trap in [^8| and find with uj± = 12n ■ 10 5 s 1 
for 53 Cr: 



^a,b,c 



0.2 

< 



hue 

Provided F is small, |Vb| can become of the order of the unperturbed level splitting tuj£, 
which is required for interactions to be relevant. 
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5. Discussion and Summary 

The prospect of realizing a quasi one-dimensional gas of interacting fermions in a trap 
will allow a comparison of experiments with emerging theories PU| , based on Luttinger 
liquid theory. These theories require a knowledge of the coupling constants describing 
the interaction. 

We have shown in detail, how a physical pair interaction between fermions in three 
spatial dimension determines three effective coupling functions V a , Vb, and V c for the 
interaction of the fermions when they are confined to a quasi one- dimensional harmonic 
trap. 

The coupling functions V a and Vb describe the forward scattering, while V c quantifies the 
backward scattering. Their respective values depend significantly on the range d^n of the 
original pair interaction, the reference length being the inverse of the one-dimensional 
Fermi wave number. 

For 2kpdsD <C 1 and distinguishable fermions, i.e., for the inter-component interaction 
in a quantum gas mixture, e.g. for the van der Waals interaction, forward and backward 
scattering contribute about equally. The effective one-dimensional potential is also short 
ranged and the equivalent of the s-wave approximation is applicable. 
For dipole-dipole interaction where d 3 n is not defined, the effective one-dimensional 
potential acquires a range din inside the trap. 

For distinguishable fermions, backward scattering becomes suppressed in comparison 
with forward scattering with increasing range of interaction. 

These results change when the fermions are identical. The contact potential ("s-wave 
scattering") then does not contribute. This has a surprising effect on the effective 
coupling functions: To a good approximation, only backscattering survives, i.e., V a and 
H, can be neglected in comparison with V c even if this backscattering coefficient is 
significantly reduced in the case of short range interactions. 

For the dipole-dipole interaction, however, it is of the same order as for the full dipole- 
dipole interaction. In reality, only the electric dipole-dipole interaction between identical 



fermionic molecules with a permanent electric dipole moment |37| can become relevant. 
Finally we mention that dipolar fermionic molecules 40 K- 87 Rb were discussed [19| as 
possible candidates for fermionic superfluidity. 
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Figure captions 



Figure 1. Reduced interaction coefficient V (s) (1220, 1214; 1200, n)/(a{V)) for "s- 
wave" scattering in the quasi one-dimensional harmonic trap. Three contributions 
at n = 1194 (Vb), n — 1206 (V c ), and n — 1234 (14) dominate. Coefficients with an 
odd sum of arguments vanish. 



Figure 2. Reduced interaction coefficients V (p) (1220, 1214; 1200, n)/(a 3 ±(z 2 V)) 
for "p-wave" scattering versus n for the quasi one-dimensional harmonic trap. 
Backscattering coefficient v} p ^ < at n — 1206 dominates. Coefficients with an 
odd sum of arguments vanish. Background coefficients as seen in figure 1 almost 
disappeared. 



Figure 3. Reduced interaction coefficients V(107, 103; 105, n)/\ Vi e ffDD(0)| for the full 
dipole-dipole interaction between atoms in different hyperfine states versus number n 
of oscillator state. Dipoles are oriented parallel to the trap axis and a filling factor of 
~ 1 is assumed. Backscattering coefficient V c at n = 109 is smaller than the forward 
scattering coefficients V a and Vb, which also dominate the other interaction coefficients. 
Coefficients with an odd sum of arguments vanish. 



Figure 4. Reduced interaction coefficients v"( non ~ s ) (107, 103; 105, n)/\ V us dd(0)| for 
the dipole-dipole interaction between identical atoms versus number n of oscillator 
state. Dipoles arc oriented parallel to the trap axis and a filling factor of « 1 is 
assumed. Backscattering coefficient vj aon at n = 109 dominates all other interaction 
coefficients including the forward scattering coefficients. Moreover, background 
coefficients are much smaller in comparison with the main peak than in figure 3. The 
inset shows the real space version of the effective one-dimensional potential Vi e ffDD 
versus z. Coefficients with an odd sum of arguments vanish. 



Appendix: Model interaction in the box 

Equations (20) and (22) allow easy numerical evaluation, but are not transparent. We, 
therefore, study a solvable model for the interaction coefficients. We consider the one- 
dimensional potential (cf. also [42] for a related analysis) 



for |z| < d and V(z) = for |^| > d, 



(A.1) 



15 

normalized according to 

/oo 
dzV(z)=g. 
-oo 

The fermions are confined to a box of length L ^> d with infinitely high walls. The 
single particle eigenstates are then 

/ 2 muz 

Wm\z) = y -y sm — —, m = 1,2,3,.... (A. 2) 

The interaction coefficients according to equation (9) are given (up to a small correction 
of order d) by 

V hox {m, p; q, n; d) = —— / dzi sm — — sm — — / dz 2 sm — — sm — — (A. 3) 
dL A Jo L L Jzi-d L L 

1 /sin(p — ri)ird/L sm(p + n)ird/L \ 

= 2L [ (p-n)nd/L ( dm+ ™ +n + 6 ^, P+q ) + ( p + n)l ,d/L dm+ ™ +n ) ' 

All coefficients, except those with 

m±p±q±n = 0, (A. 4) 

vanish due to local translational invariance inside the box. But near the Fermi surface 
(m,p, q, n ~ N ^> 1) only the three processes shown in equation (A. 3) remain: two with 
small and one with large momentum transfer Ap as was discussed in Section 2. 
In the notation of equation (11), this can be written as 

+ V b 6 + V C 8 (A.5) 

Again, for processes near the Fermi edge and d<L we approximate 

AAfezg"), (A,) 




sin(2A%d/L)\ 
2Nird/L J ' 



(A.7) 



The coefficients depend on the differences between the indices only (V^b) or are 
constant with respect to the indices (V c ). We note in passing that the interaction 
coefficients Vb ox (m,p; q,n; d) show perfectly the behaviour which allows application of 
the Bosonization method, see JftJ . 
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In the limit d — > the potential ( |A.1| ) becomes a contact potential (^(z) with pure 



s-wave scattering and all three types of coefficients are equal and constant 

K,b, c - K ( t = ^ (A-8) 

However, when the range d of the potential ( |A.1| ) increases and becomes larger than the 
average distance L/A of the atoms in the box the backscattering coefficient V c decreases 
while V^b remain constant as long as d <C L. 

As discussed in Section 3 above, any fully symmetric part common to all interaction 
coefficients does not contribute to the interaction V3 according to equation (|TD and 
may therefore be removed. Doing so we subtract the fully symmetric potential 
^box(^,p; q, n; d = 0) (|A.8|) from Vb ox (m,p; q, n; d) (|A.5| , |A.6|JA.7|) with the result that 



only a backscattering term which we denote by remains 

V hox (m,p; q, n; d) -> V c {p) 5 m+qjP+n (A. 9) 

S m(2nNd/L) \ 

2nNd/L - 1 5 - +q , P+n (A.10) 




~ -^Lkld 2 5, 



3L 



m+q,p+m 



where the latter approximation is valid for k-pd 1 (kp = ttN/L). The sign of K (p) 
is opposite to the sign of in accordance with the results for the harmonic trap, 
compare figures 1 and 2. 

In equations (|2T| ) and (^3) we found in the case of the harmonic trap for the scaling 
of the coefficients with the particle number A^: V^ c ~ A _a3 and ~ A 5 . To 
compare with equations ( |A.8| ) and ( |A.9| ) where the coefficients are given in terms of L 
and k-p we recall that in the harmonic trap kp and L are proportional to y/2N. We 
substitute these relations into the box coefficients and get 

K ( ,b,c ~ N- 1/2 , K (P) ~ N 1 ' 2 . (A.ll) 

We conclude that in spite of the crude character of our model the agreement of the 
cofficients in the box and in the harmonic trap is quite good. The background 

terms for the momentum non-conserving processes in the harmonic trap are, of course, 
absent in the box. 
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